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Tsoukias, NlkoJaos M, and Steven C. George. A two- 
compartment model of pulmonary nitric oxide exchange 
dynamics. J. Appl. Physiol. 85(2):653-666, 1998.—The rela¬ 
tively recent detection of nitric oxide (NO) in the exhaled 
breath has prompted a great deal of experimentation in an 
effort to understand the pulmonary exchange dynamics. 
There has been very little progress in theoretical studies to 
assist in tire interpretation of the experimental results. We 
have developed a two-compartment model of the lungs in an 
effort to explain several fundamental experimental observa¬ 
tions. The model consists of a nonexpansile compartment 
representing the conducting airways and an expansile com¬ 
partment representing the alveolar region of the lungs. Each 
compartment is surrounded by a layer of tissue that is 
capable of producing and consuming NO. Beyond the tissue 
barrier in each compartment is a layer of blood representing 
the bronchial circulation or the pulmonary circulation, which 
are both considered an infinite sink for NO. All parameters 
were estimated from data in the literature, including the 
production rates of NO in the tissue layers, which were 
estimated from experimental plots of the elimination rate of 
NO at end exhalation (E N0 ) vs. the exhalation flow rate (V E ). 
The model is able to simulate the shape of the NO exhalation 
profile and to successfully simulate the following experimen¬ 
tal features of endogenous NO exchange: I) an inverse 
relationship between exhaled NO concentration and Ve,. 2\ 
the dynamic relationship between the phase III slope and Ve, 
and 3) the positive relationship between E^o and Ve. The 
model predicts that these relationships can be explained by 
significant contributions of NO in the exhaled breath from the 
nonexpansile airways and the expansile alveoli. In addition, 
the model predicts that the relationship between E N0 and V E 
can be used as an index of the relative contributions of the 
airways and the alveoli to exhaled NO. 

mathematical model; elimination rate; phase III: single exha¬ 
lation 


THE IMPORTANCE OF NITRIC Oxide (NO) in the physiology 
of lung function has steadily increased in the past 
decade. Experimental work has dominated this influx 
of information and has provided critical Information 
related to the cellular source of NO in the lungs, the 
specific physiological functions, and the detection of 
endogenous NO in the exhaled breath. On the basis of 
the ability of NO to relax smooth muscle, there has 
been interest in using exogenous NO delivered by 
inhalation as a therapy in such diseases as pulmonary 
hypertension, acute respiratory distress syndrome, and 
bronchial asthma. In addition, there has been interest 
in using endogenously exhaled NO levels as a noninva- 
sive index of pulmonary inflammation (13). This inter¬ 
est has generated a need to understand the exchange 
dynamics of NO within the lungs. 
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The exchange dynamics of NO in the lungs are different 
from those of other previously studied endogenous gases 
such as O 2 and C0 2 because of specific differences in 
physical and biochemical properties. I) NO, with 11 va¬ 
lence electrons, is a free radical and is reactive with a 
number of different endogenous substrates, including 0 2 , 
superoxide, thiols, and metalloproteins. The presence of 
this array of substrates results in a half-life ( 4 .) of free NO 
in vivo that is ~0.5-15 s, depending on the specific 
physiological state (5, 16). 2) NO is actively produced 
basally and in response to various stimuli, such as cyto¬ 
kines, by the ceils that comprise the alveolar and airway 
regions of the lungs. 3) NO binds avidly to hemoglobin in 
the blood, the kinetics of which are much different from 
those of CO 2 and O 2 , but may play an important role in the 
local control of oxygenation (12). Theoretical modeling has 
historically played an important role in our understanding 
of pulmonary gas exchange dynamics; however, because of 
the unique biochemical features of NO, new models must 
be developed to complete our understanding of NO ex¬ 
change. The lone attempt to model NO exchange in the 
lungs lumped the entire lung into a single compartment 
that effectively represented the alveolar region (II). This 
model was a starting point for modeling NO exchange 
dynamics and provided some initial insight into NO ex¬ 
change. However, much further work is needed, inasmuch 
as a single-compartment model is not capable of describing 
the growing experimental evidence that exhaled endog¬ 
enous NO is derived from an alveolar and an airway 
source. The alveoli and airways are two distinct regions of 
the lungs with many different features. From a gas- 
exchange perspective, perhaps the most important is the 
fact that the alveolar volume is expansile and the airway 
volume is relatively nonexpansile. The goal of this study is 
to develop a relatively simple two-compartment model of 
NO exchange dynamics that can /) describe the basic 
features of NO exchange dynamics previously ohserved 
experimentally by us (26) and by others (20, 24), 2) provide 
initial insight into the exchange properties of exogenous 
NO, and 3) provide direction for future experimentation. 

MODEL DEVELOPMENT 

General Structure 

A model of the human lungs was developed (Fig. 1 ) 
that includes enough detail to describe the following 
basic features of NO exchange observed experimen¬ 
tally: 1) breath holding before exhalation creates an 
initial spike in NO concentration during exhalation of 
the airway volume, Z) exhaled NO concentration is an 
inverse function of exhalation flow rate (Ve), and 3) 
elimination rate of NO from the lungs (E NQ ) is a positive 
function of V E . These basic observations are consistent 
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TWO-COMPARTMENT MODEL OF NO EXCHANGE 



Fig. 1. Schematic of 2-compartment model for nitric oxide (NO) 
pulmonary exchange. First compartment represents relatively nonex- 
pansile conducting airways; second compartment represents expans¬ 
ile alveoli. Each compartment is adjacent to a Layer of tissue that is 
capable of producing and consuming NO. Exterior to tissue is a layer 
of blood that represents bronchial or.pulmonary circulation and 
serves as an Infinite sink for NO, Ve and V), expiratory and 
inspiratory flow, respectively; Ce and Ci, expiratory and inspiratory 
concentration, respectively; C.nr and C&iv. airway and alveolar concen¬ 
tration, respectively; V alr and V aiv , airway and alveolar volume, 
respectively; -/ t;g ,ai r and -/ l:&a | v , total flux of NO from tissue to air and 
from alveolar tissue, respectively; t, time; V, volume. 


with the alveolar and the airway regions of the lungs as 
sources of exhaled NO. Hence, the model (Fig. 1) 
consists of two main compartments; 1) a rigid or 
nonexpansile compartment representing the conduct¬ 
ing zone or airways [trachea through generation 17 as 
defined by Weibel (27)] and 2) a flexible or expansile 
compartment representing the respiratory bronchioles 
and alveolar region (generation 18 and beyond). Both 
compartments are surrounded by a layer of tissue 
representing the bronchia] mucosa in the airway com¬ 
partment and the alveolar membrane in the alveolar 
compartment. Blood representing the bronchial circula¬ 
tion and the pulmonary circulation is distal to the 
tissue in the airway compartment and alveolar compart¬ 
ment, respectively. Cells present in the airway mucosa 
and the alveolar membrane are capable of producing 
endogenous NO (8, 23); hence, we assume that NO is 
produced at a constant rate per unit volume of tissue by 
the tissue surrounding the airway and alveolar compart¬ 
ments. The endogenously produced NO can follow one 
of three paths: 1) consumption through reaction with 
substrates within the tissue compartments, 2) diffusion 
toward the pulmonary or bronchial circulations, where 
it reacts instantaneously and irreversibly with the 
hemoglobin of the blood, or 3) diffusion toward the 
airstream, where it evaporates and enters the alveolar 
volume or airway volume. Consequently, there will be a 
net flux of NO between the tissue and the airstream in 
both compartments. The direction of the flux depends 
on the relative concentrations of NO in the tissue and 


gas phases. Hence, NO concentration in the exhaled 
breath will depend on two additive mechanisms: 1) the 
exchange of NO in the alveolar compartment and 2) the 
conditioning of the alveolar gas as It is convected 
through the airway compartment. On the basis of this 
general structure, a series of mass balances on the 
tissue and gas phase of the airway compartment and 
the alveolar compartment produce the governing equa¬ 
tions for the model. 

Airway Compartment 

Tissue phase. The airway compartment is modeled as 
a cylindrical tube, and the surrounding tissue (Fig. 2) is 
modeled as a homogenous compartment of uniform 
thickness (/.. „ lr ) with aqueous physical properties. NO 
is produced uniformly in position and constant in time 
by the tissue at a rate S tialr (mol ■s"' 1 -cm" 3 ). NO reacts 
with several substrates including O 2 , superoxide, thi¬ 
ols. and metalloproteins such as guanyiate cyclase (5); 
hence, NO is consumed by chemical reaction at a rate 
R ta i r (mol-s" 1 -cm -3 ). The reaction with 0 2 is second 
order in NO concentration, and this reaction is the 
rate-limiting step in the reaction with endogenous 
thiols such as glutathione (14). At physiological NO 
concentrations, only the reactions with superoxide and 
metalloproteins are significant (2); both are first order 
in NO concentration. Hence, the R t alr can be expressed 
as x, z), where k is an overall first-order rate 

Constant (s“ J ), C t (f, x, z) is the concentration of NO in the 
airway tissue, x is the distance from the bronchial blood 
compartment, zis the axial position, and tis time (s). 



Fig. 2. Schematic of tissue layer in airway compartment. NO is 
produced from cellular sources at a constant rate in position and time 
(St,air)- NO can then be transported by diffusion and enter bronchial 
circulation or gas phase of airway lumen, or it can be consumed by 
chemical reaction that is 1st order in NO concentration. Solid line, 
approximate radial concentration profile of NO in blood (b), tissue (t), 
and gas (g) phases under ambient inspired conditions. Radial gradi¬ 
ent in NO concentration at tissue-air interface is proportional to flux 
of NO into airstream. A,air. thickness of tissue surrounding airway 
compartment; partition coefficient of NO between tissue and air; 

Cw, concentration at airway-tissue interface; C t , concentration in 
tissue. 
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The transport of NO from the site of production to the 
bronchial blood or the airstream occurs by molecular 
diffusion described by Fick’s first law. The axial and 
angular gradients of NO concentration in the tissue are 
considered negligible. Hence, the transport of NO in the 
tissue can be sufficiently described by a one-dimen¬ 
sional diffusion equation. The small tissue thickness in 
comparison with the airway radius allows the use of 
simple Cartesian coordinates. Because of the fast reac¬ 
tion of the NO with hemoglobin present in abundance 


Once the concentration profile in the tissue is known, 
using Fick's first law of diffusion, one can easily deter¬ 
mine that the flux of NO from the tissue to the air 
(-fi: MlS ]r) is a linear function of C w 

ac t i 

-4:g,air = — A ~ a ~ b * C w (4a) 

-k^h.air 

where 


a = — 


J It,air/3 

/St.airj 

It,air/P 

UJ 


, 

• air/3 -|- q — 


St.a 




m 


in the blood, the concentration of free NO at the 
interface between the blood and the tissue will be very 
close to zero. For the interface between the airway 
lumen and the tissue, we assume instant thermody¬ 
namic equilibrium governed by Henry’s law. Under 
these assumptions, a differential mass balance of NO in 
the tissue produces the following second-order partial 
differential equation 


a 2 C t 
1? + St ’ air 



( 1 ) 


with the following boundary conditions: C t (q 0) = 0, C t (f, 
Aair) = Cw(f, z). D t is the molecular diffusion coefficient of 
NO in tissue (3.3 x ICO 5 cm 2 /s) (16), and C w is the concen¬ 
tration of NO at the interface between the airway and 
tissue. C w (t, 2 ) will depend on the concentration in the gas 
phase; as a result, Eq. 1 must be solved simultaneously 
with tile NO mass balance in the gas phase of the airway 
compartment (see below). The solution to Eq. 1 can be 
simplified by assuming that C t is in the steady state. This 
approximation is valid if tire times of inspiration and 
expiration are much greater than ~0.6 s (the time for tire 
concentration profile in the tissue to reach 90% of its 
steady-state value). This is due to the relatively small 
thickness of tissue relative to the rate of diffusion, and a 
proof is provided in appendix a. With the assumption of 
steady state in the tissue, Eq. 1 reduces to 

d 2 Ct 

A-^r + S tialr -AC t = 0 (2) 

The solution of Eq. 2is 


b 


A 

P 


l 

air/p — i-t.alr/P 


-|- g»“^trafr/P) 


(4c) 


Because of the very low solubility of NO in water and 
tissue, the radial transport of NO is not limited by 
diffusion in the gas phase (appendix b) (6). Hence, Eq. 4 
becomes 


rfcg.air — a b * (k t:s C a | r ) (5) 

where A l:g is the partition coefficient of NO between the 
tissue and the air at 37°C and C air is the bulk gas 
concentration of NO in the airway lumen. From Eq. 5, 
the flux of NO per unit airway surface area and per unit 
time between the airway tissue and the airway lumen 
is a linear function of the bulk gas concentration. As the 
concentration in the bulk gas of the airway lumen 
increases, the amount of NO that is consumed by the 
pulmonary blood or through reactions with substrates 
in the airway tissue increases; thus galr decreases. 
For Cai, greater than a/(b\ t: g), the consumption over¬ 
comes the endogenous production rates, resulting in 
net flux of NO from the airway to the tissue. 

Gas phase. The airway compartment is modeled as a 
cylinder of constant total volume Vaj r and total surface 
area A s . alr . As described above, the absorption of NO by 
the airstream per unit surface area at any given point 
of the airway is described by Eq. 5. If we assume plug 
flow (i.e„ no radial gradient in the velocity or concentra¬ 
tion) , then a differential mass balance of NO over a unit 
volume will give 


. dC„ 


dz 


T A:g,alr 




dt 


(6a) 


C t (x) = «(e tf P - e-**) + (1 - 


. . which can be transformed to a more convenient form by 

substituting dV = dzA c 


where 


git.air/P 


c w + 


-WlS w 


-1 + a-WP js xMr \ 

£^t,alr/P — £^£t,alr/p \ A j 


(3b) 


P - 



(3c) 


-V- 


ac. 


air 


3V 


+ T —- 

Cg,alr a 


aCpir 

dt ' 


(6b) 


where V is volumetric flow rate of air (Vj for inspiration 
and V E for expiration), A s is the wall surface area per 
unit axial distance, A c is airway volume per unit axial 
distance, and V is axial position in units of cumulative 
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volume. The ratio AJA K (surface area per unit airway 
volume) increases in the human lung as the radius of 
each airway branch (r) decreases with generation num¬ 
ber {\iA v — 2/r). Hence, a more realistic description of 
the airway compartment is not as a perfect cylinder 
with constant A s /A c but, rather, as a cylinder whereby 
AJA C increases with increasing axial position (or de¬ 
creasing radius). To preserve an analytic solution to the 
model equations, A s /A c Is approximated by a linear 
function of V (7? = 0.94) with use of the airway dimen¬ 
sion data of Weibel (28): AJA Q = aj + by for inspiration 
and AJA q = % + by for expiration (a E = a x + by^, 
b£ — — b[). Hence, the total surface area of the airway 
compartment (A . 3lr ) is described by 



Equation 6 is subject to the following arbitrary initial 
and boundary conditions: C 3lr (0, V) = /(V) and Cau-fc 0) = 
git), where /(V) is the concentration profile of NO in the 
upper compartment before the inspiration (V — Oat the 
mouth) or expiration (V = 0 at the entrance of the 
alveolar compartment) and g(t) is the inspired concen¬ 
tration or the alveolar concentration for inspiration (I) 
or expiration (E), respectively. The solution to Eq. 6 foe 
constant V is 


phase ( Jt;g.aiv) is a linear function of the bulk gas phase 
concentration (see Eq. 5). Then the net total flux (mol/s) 
of NO from the alveolar tissue (-/fg.aJv ^s.uiv, where ^4s ia iv 
is the surface area of the alveolar region) can be 
described by a combination of two terms 

“^t:g,alv ■'^s.alv ^app,alv Hl^qC 3 i v (9) 

where S appa iv (mol/s) is the apparent production rate of 
NO is the tissue (defined as the flux of NO into the 
alveolar compartment if the concentration of NO was 
zero in the compartment) and Dl^q is the diffusing 
capacity (mo! 'S' 1 •mol" 1 -cm 3 ) of NO in the alveolar 
region. The method for determining the alveolar flux 
contrasts with that used in the airway compartment 
only in the manner of determining the coefficients of 
the linear function (i.e., S appialv and Dl no in Eq. 9). In 
the alveolar compartment the coefficients can be deter¬ 
mined experimentally rather than calculated from the 
geometric characteristics of the alveolar tissue and the 
endogenous production rate, as was done in the airway 
compartment. Dl n0 has been experimentally measured 
(3, 7) and is equal to ~2,10Q mol-s" 1 -mol -1 -cm 3 , and 
S a pp,aiv can be determined from the steady-state alveo¬ 
lar concentration (see Eq. 12). 

A differential mass balance for NO in the alveolar 
compartment (valid for inspiration and expiration) is then 


Caj r (r, V) = 


/(V - Vt) - 


b\ t . f 


'“'**•»* -I- —— (8a) 

0\.g 


for case I (t < V/V) and 


Carte V) = 


Lgte — v/v) 


bK, 


-gJ 


a(v/V) + __ (#/j) 

OVg 


for case II (t > V/V), where 


y(f, V, j) =aj+ b/S - Vt/2) (Sc) 

S(V, j = aj + bjVI 2 J = 1, E (8d) 

Case I descri bes the convective emptying of the airways 
at the beginning of inspiration or expiration; case II 
describes the convection of the alveolar gas or inspired 
gas through the airway compartment during expiration 
or inspiration, respectively. For the linearly increasing 
and decreasing flow rate maneuvers, V E is not constant. 
However, an analytic solution is still attainable and is 
presented in appendix c. 


Vaiv(i) 


dC aJv 

"dT 


UO) 


('Npp.alv HLjsj^jC-j v ) -p Vj(Cg| rien d C^y) 

subject to the following initial condition: Caiv(O) — Calv,0 • 
VaivW is the volume of the alveolar compartment, and 
Cair.and is the concentration of NO in the air that exits 
the airway compartment. The volume change due to 
exchange of the respiratory gases has been neglected. 
Equation 10 is valid for inspiration and expiration 
(Vi = 0) as well as for nonconstant flow rates. 

The analytic solution of Eq. 10 for constant flow rates 
and for constant Cau^d is as follows 

IVaivPF - 1 

cwo - c+m [y-^j 

Sapp.alv ~f VlCaif.etiti 

DL N0 4- Vj 


alv(9) 


K+r 


V alv (/) 


(11a) 


Alveolar Compartment 

Gas phase. In modeling the alveolar region (Fig. 1), 
we assume a well-mixed compartment of variable vol¬ 
ume IVajvMJ- The NO concentration in the alveolar gas 
[C alv (t)] is uniform in position but is time dependent. 
NO can enter or leave the compartment through convec¬ 
tive flow during inspiration or expiration, respectively, 
and can exchange with the alveolar tissue by diffusion. 
After an analysis that is identical to that presented for 
the airway tissue compartment, it can be shown that 
the flux of NO between the alveolar gas and the tissue 


for case I (inspiration) 


1f Pilv(O) 


V alv (t) 


ValvP 


+ 


^app.NO 


Dt, 


NO 




V^O) 


(lib) 


for case II (expiration), where ^ = Dl n0 Nj (J = I or E), 
and 

C alv (i) = C a j v (0) e'l DL NoV alv «k 


^app.alv 


Dl, 


NO 


[1 


-iDLKoflljvMJtj 


(lie) 
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for case III (breath hold). Solution of Eq. 10 for linearly 
time-dependent expiratory flow rate is derived in appen¬ 
dix d. 


Parameter Estimation 

Table 1 summarizes the values for key parameters 
used in the model simulations. The A t:g is calculated 
from Henry’s law constant for NO between air and 
water (25.8 X 10 6 mmHgNO-mol NO -1 -mol H 2 0) (18a). 
The reaction constant k is calculated from the f w of NO 
in tissue. For first-order reactions, k = In (2)/t Y „ The 
values of L ri that have been reported in the literature 
range from 0.5 to 15 s, depending on the tissue and local 
environment (15, 16, 28). For our simulation, we chose 
an intermediate value of 4 s. The volume and the 
surface area of the upper compartment are calculated 
from Weibel’s anatomic data (28) between the trachea 
and generation 17. The coefficients of the linear regres¬ 
sion of the ratio were chosen to satisfy Eq. 7. 

S ap p.aiv and S r ae are more challenging to estimate. 
Sapp.aiv can be estimated from the diffusing capacity of 
NO and the expired steady-state alveolar concentration 
(C a t,ss, Eq. 10) if Ca[ V . ss is known 

^app.alv ^1-NO^alv.ss ( 1 ^) 

Hyde et al. (11) used the end-exhaled concentration of 
NO after a 15-s breath hold as an estimate of C bJv . ss . 
However, this approximation is valid only if the contri¬ 
bution from the airways is neglected [recall that the 
model of Hyde et al. (11) had only an expansile alveolar 
compartment]. An alternative technique, which ac¬ 
counts for the contribution from the airways, is to use 
E N o from the lungs, as described in the companion 


Table 1. Parameter values for two-compartment model 


Parameter 

Value 

Units 


Tissue 

0.0416 

ranffmol 

»c 

3.3E-05 

cm7s 

k 

0.173 

s -1 


0.01 

cm 

St.air 

0.55 

pmol*s -1 ’cm -3 

^s.alr 

Airway 

5,100 

cm 2 

v air 

200 

cm 3 

ai 

12 

cm -1 

bl 

0,333 

cm 2 

Vaiv 

Alveolar 

2,300-5,800 

cm 3 

D4io 

2,100 

mats -1 -mol -1 -cm 3 

Sapp.aiv 

670 

pmal/s 


\t; gl Partition coefficient of NO between tissue and air; D,. molecu¬ 
lar diffusion coefficient of NO in tissue: k. overall 1st order rate 
constant; A,. a i r , thickness of tissue surrounding airway compartment: 
St.air. rate of NO production by tissue in airway compartment; A Si3Jr . 
surface area of airway compartment: V alr , volume of airway compart¬ 
ment; ai and A, constants for inspiration; V„i v , alveolar volume; 
Dlno, diffusing capacity of NO in alveolar region; S jp p. a ] U , apparent 
production rate of NO in alveolar tissue. 


article (27) 

— Vli.tcCalv.cG fft:g,alr Tls air (13a) 

3: g . a lr = 77—7,-g.alrdV 

v air u 

where ,7 cgair is the average (over axial position) J t:gaJr 
and V Eiee is V E at end exhalation. Hence, an estimate 
of C a j„ iSS .(Cgiyeg) is attained from the slope of a plot of 
E n o vs. V E ee . We used the elimination rate data from a 
representative subject ( subject 6) in the companion ar¬ 
ticle (26) for the model simulations. For subject G.C^v.ss 
was equal to 8.11 parts/billion (ppb), and thus S app , 3 i„ 
was estimated as 6.7 X 10" 10 mol/s. 

au can also be estimated from Eq. 13, with the 
intercept of the E N0 vs. V Eit:e plot as an estimate of 
Tog,airmail- St,air is then chosen using an iterative 
scheme to satisfy Eqs. 4 and 13 simultaneously, for 
which St.air is the only unknown parameter. In the case 
of subject 6, -/uftairA.air was equal.to 2.33 X 10 -!1 mol/s, 
which produced a value for S ta i r of 5.5 x 10 ~ 13 
mol-s - 1 -cm -3 . 

The model equations can be solved for different 
breathing patterns of constant Ve as well as Vg that 
change linearly in time (see appendices c and d). The 
initial alveolar volume, inspired and expired volumes, 
and inspired and expired flow rates are specified by the 
user. Unless otherwise noted, inspiratory conditions for 
all simulations were from functional residual capacity 
to total lung capacity (alveolar volume 2,300 and 5,800 
ml, respectively). Functional residual capacity and 
total lung capacity were estimated from the vital 
capacity (VC) of subject 6 (9). The control inspiratory 
conditions also included the following parameter val¬ 
ues: Vi = 2,000 cm 3 /s, C a i v0 = 8.3 ppb, and Cj = 15 ppb, 
where C^q is initial alveolar concentration and Q is 
inspiratory concentration. C £jvCl can be estimated by 
first setting its value to zero and then simulating 
several rebreathing maneuvers (tidal volume = 500 
ml, respiratory rate = 12 breaths/min) by setting the 
initial concentration to the end-expiratory concentra¬ 
tion of the previous maneuver. A steady-state value is 
quickly achieved, which can serve as C.^.q for the 
remaining simulations. 

RESULTS 

Simulation of the Single Exhalation 

Figure 3 depicts four different experimental single¬ 
exhalation maneuvers of the representative subject 
(26) and the model prediction of these maneuvers with 
the same set of parameters (Table I) in each case. 
Recall that no particular optimization procedure has 
been performed to simulate the experimental profiles. 
All parameters were estimated from the literature or 
from the elimination rate data from this particular 
subject. In Figure 3, A and B, the subject exhales at a 
constant low flow rate (control) and a constant high 
flow rate, respectively. In Fig. 3 C, the control maneuver 
is performed after a 15-s breath hold, and in Fig. 3D V E 


PM3006723688 


Source: https://www.industrydocuments.ucsf.edu/docs/qgbj0001 



658 


TWO-COMPARTMENT MODEL OF NO EXCHANGE 


Fig. 3. Model simulation (solid lines) 
of 4 different experimental exhalation 
NO profiles in which Ve is indepen¬ 
dent variable. A: control maneuver, 
which represents a constant Ve of—200 
cmVs. &. high flow rate maneuver in 
which Ve is still constant but is in¬ 
creased to —450 cra 3 /s. C: preexpira- 
toiy 15-s^breath hold followed by a 
constant Ve. of —200 cm%. D\ linear 
decrease of Vg in time during course of 
exhalation with a slope of approxi¬ 
mately —40 cm 3 -s -I *s and an average 
Ve of —300 ciriVs. Ce, expiratory con¬ 
centration; ppb, parts/billion. 


A 

so - 



Time (sec) 

B 


20 T 



C 




Time (sec) 


is decreasing linearly with time on the basis of the flow 
signal from this particular profile. 

At a constant V E (Fig. 3. A and B) the model 
replicates the initial rise in exhaled NO concentration 
(beginning at Cj as the airway compartment is exhaled 
(first 200 ml or ~1 s). This rise represents an axial NO 
gradient due to the absorption of NO by the airstream 
from the airway tissue. After this initial peak, there is a 
step change as air from the alveolar compartment 
reaches the exhalate. This value (10 ppb) is larger than 
the steady-state alveolar concentration (8.11 ppb) and 
the initial or preexpiratory concentration (8.3 ppb) due 
to absorption of NO from the airway tissue as the 
airstream is converted through the airways. The model- 
predicted and experimental expiratory concentration 
(C E ) then steadily declines over the remaining portion 
of the exhalation. After a 15-s breath hold (Fig. 3 C] the 
model predicts a large initial peak due to accumulation 
of NO in the airway compartment followed by a similar 
decline, although less pronounced, in C E over the 
remaining exhalation. When V E decreases linearly in 
time (Fig. 3D) the model is able to simulate the initial 
fall in C E during phase III followed by a steadily 
increasing rise. 

Flow Rate Dependence of the Average- 
Phase III Concentration 

In Fig. 4, average phase III concentrations with 
respect to time [C E (f)] or volume [C E (V)| are plotted 
as a function of V E over a wide range of constant V E . 
C E (t) is estimated over a specified time interval (exhala¬ 
tion time of 2-8 s) and C E (V) over a specified 
volume interval (exhalation volume of 20-60% 
of VC), as described in the companion article (26). For 


very low or very high V E , where phase III does not 
include the above intervals, the maximum allow¬ 
able interval is used. The experimental data from 
the companion article are also reported for reference. 
There is a good agreement between our experimental 
data and our simulations for intermediate flow 
rates. Over the entire flow rate range there is an 
exponential decrease for C E (f) and C E (V) with V E . This 
result is consistent with the experimental data (end- 
exhaled NO concentration) presented by Silkoff et al. 
(24) over the range 4.2-1,500 cm 3 /s. For low V E (<500 
cm 3 /s), C E (f) has a slightly larger value. In the high-V E 
region (>500 cm 3 /s), C E (V) is slightly higher and has a 
small positive slope. In the intermediate region the 
slope of C^t) vs. flow (—0.004 ppb-mF 1 -s -1 ) is steeper 
than for C E (V) (—0.002 ppb-ml -1 -s _l ), which is in 
agreement with our experimental data (see Fig. 4 in 
Ref. 26). 
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Fig. 4. Dependence of average NO concentration in phase III (C E ) 
with respect to time ( t ) and volume (V) over a wide range of V E , 
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Phase III Slope 

Figure 5 presents data from the simulation of lin¬ 
early changing flow rate profiles. The normalized phase 
III slope (slope in phase III normalized by the average 
concentration in phase III) in time [SnuyoWL is plotted 
as a function of the normalized slope of Ve [Sm,v E (f): see 
Fig. 6 of Ref. 26], To generate the data, the model 
simulated the exhalation maneuvers of subject 6, match¬ 
ing the experimental flow rate profile of each maneu¬ 
ver. The simulations are performed for the breath-hold 
and the non-breath-hold maneuvers. The experimental 
data for subject 6 are also shown for comparison. The 
model is able to simulate the inverse relationship, 
including the negative intercept at a constant V E . The 
inverse relationship is due to the nonnegligible contri¬ 
bution from the airways to exhaled NO and, hence, the 
flow rate, dependence. As V E decreases, C E Increases. 
Thus, if Ve decreases in time, Ce will increase in time 
and affect the slope of phase III. The slight negative 
slope at constant V E is generated from continuous gas 
exchange in the alveolar region. In addition, in a 
fashion similar to the experimental data, the intercept 
increases after a 15-s breath hold. This phenomenon is 
due to a decrease in the alveolar concentration toward 
the steady-state value during the 15-s breath hold (see 
Inspiratory Conditions). 

Elimination Rate 

In Fig. 6 we present E^q at the end of expiration, 
estimated as the product of end-expiratory exhaled 
concentration and Vg. The data are derived from model 
simulations of constant Ve (range 50-1,000 cm 3 /s) 
exhalation maneuvers with or without a 15-s breath 
hold. For V E > 500 ml/s only data from the breath-hold 
experiments were used to ensure an elapsed time of 
5:8 s from end inspiration (see Inspiratory Conditions). 
Least squares linear regression over the resulting data 
reveals a highly linear {R i = 0.99) dependence of E N0 
with V E . The slope and intercept of the linear regres- 



Fig. 5. Experimental and model-predicted relationship between nor¬ 
malized phase III slope of exhaled NO concentration in time 
and exhalation flow rate JA n -.i, b (/J|wiih and without a 15-s 
breath hold. 



(ml/s) 


Fig. 6. Model-predicted relationship between elimination rate of NO 
at end exhalation (Exo) and Ve. E»o is defined as product of 
end-exhaled concentration and Ve. Positive slope is consistent with a 
nonzero alveolar concentration of NO; nonzero intercept is consistent 
with a nonnegligible flux of NO from airway tissue. 


sion line are 8.12 ppb and 0.589 X 10 -6 ml/s, respec¬ 
tively, which are in close agreement with the experimen¬ 
tal values from subject 6 (8.11 ppb and 0.592 X 10~ 6 
ml/s). The nonzero slope and intercept are due to 
significant contributions of exhaled NO from the expan¬ 
sile alveoli and the nonexpansile airways (Eg. 13). 

Inspiratory Conditions 

In Fig. 7 the model-predicted effect of the inspired 
conditions on the exhalation profile of NO is examined. 
Although this represents an extrapolation of the mode], 
this type of prediction can be used to guide future 
experimentation. Figure 7.4 presents the model's predic¬ 
tion for the NO single-exhalation profile. Three differ¬ 
ent values for the inhaled concentration are shown: 0, 
15, and 30 ppb. For Q = 0 the model predicts an initial 
rise from 0 ppb as the conducting airway space is 
emptied, then a more gradual rise to the same plateau 
value seen previously for C; = 15 ppb. The result Is a 
positive slope for phase III. In addition, the model 
predicts obliteration of the peak in phase I after 
inspiration of NO-free air. For C r = 30 ppb, there is, 
again, an initial rise in C E as the airway compartment 
is emptied, which is then followed by a steeper decline 
to the end-exhalation value compared with Ci — 15 ppb. 
Hence, the model predicts that ambient concentration 
can significantly affect the peak and the slope of phase 
III during exhalation. 

Figure IB presents the effect of Vj on the exhalation 
profile when Vj, C;, and V E are held constant at the 
control values (Q = 15 ppb, Ve = 250 cm 3 /s, V- — 3,500 
ml). The model predicts that the slope of phase III 
becomes progressively smaller as V r decreases from 
2,000 to 100 ctrfVs, yet the end-exhalation concentra¬ 
tion remains unchanged. In addition, Vi affects the 
peak in phase I; the maneuver of the slowest V : 
produces the highest peak. 

Figure 7C presents the effect of Vj on the exhalation 
profile. The model predicts that only the slope of phase 
III becomes progressively larger as Vi decreases from 
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A 




C 



Fig. 7. Effect of inspired conditions (flow rate, volume, and concentra¬ 
tion) on shape of exhalation profile and alveolar concentration. A: 
exhalation profile at 3 different Ci. B\ exhalation profile at 3 different 
Vi. C: exhalation profile at 3 different Vr. 


3,500 to 500 cm 3 , while the peak In phase I remains 
unaffected. 

The above results can be understood by examining 
the alveolar concentration of NO under these condi¬ 
tions. Figure 8 demonstrates the alveolar concentration 
of NO at end inspiration (C alVie! ) under different condi¬ 
tions of C] and Vj. C a i ViS5 is 8.11 ppb. For Q in excess of 
~5 ppb, there is an increase in C Hlvel above C alViM . 
Hence, during exhalation there is a net negative flux of 
NO {e.g., toward the pulmonary blood) in the alveolar 
compartment as DL N0 C alv > S appia!v (Eq. 9). The result is 
a progressively decreasing C a[v , and hence C E , during 
exhalation until reaches Caiv.ss- At this point,. C E 
reaches a plateau value of 10 ppb (for the control V E ). 
The opposite scenario occurs for Q less than —5 ppb. 
Under these conditions, C a [ v . e i is below C„| ViSS , and 
DL N0 Cai v < S a pp ia j y . The result is an exponential rise in 



Vj (ml/s) 

Fig. 8. Alveolar concentration at end inspiration (C^ v .ei) as a function 
of Ci and Vi. 


and Ce daring exhalation. The time constant for 
this exponential decline or rise is equal to Va[ V (£)/DL N0 
or —1-3 s over the course of exhalation. Thus the model 
predicts that an exhalation time of ~8 s will achieve an 
alveolar concentration that differs by < 1 % of the steady- 
state concentration. In addition, the difference between 
Caiv.ei and for any given inspired concentration is 
increased with Vj. As a result, phase III slope is also 
affected by Vj (Fig. 7.S). 

Concentration Profile in the Airway Tissue 

The steady-state concentration profile in the airway 
tissue, as predicted by the model, for different inhaled 
NO concentrations is shown in Fig. 9. The radial 
distance (x) is considered zero at the blood-tissue 
interface. For Q less than —418 ppb (derived from zero 
net flux or the ratio alb from Eq, 4 ), which is usually the 
case for inspiration of ambient air, there is a net 



Fig. 9. Model-predicted steady-state radial concentration profile in 
airway tissue as a Function of Ci- A: low concentrations [0 < Cj < 1 
part/million (ppm)]. B: high concentrations (Q = 50 or 100 ppm). 
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positive flux of NO or absorption of NO by the air- 
stream. For Ct greater than ~418 ppb (i.e., when 
exogenous NO is inspired for clinical application), the 
net flux of NO is negative, or there is desorption of NO 
from the airstream to the airway tissue. At these higher 
concentrations the tissue concentration profile is ap¬ 
proximately linear when 4 is 4 s. Under these circum¬ 
stances, the consumption rate of NO within the tissue 
is small compared with the flux of NO from the 
airstream. 

Exogenous NO 

In Fig. IQ A we present the model predictions for the 
fate of exogenous NO at concentrations that are repre¬ 
sentative of clinical situations [C; = 10 parts/million 
(ppm)]. A single rebreathing maneuver.of the. tidal 
volume (V] = 500 ml) is simulated, while Vj and V E are 
held constant (V t and V E = 250 ml/s). The model pre¬ 
dicts that 55.1% of inspired NO (90% of the NO that 
actually reaches the alveoli) is absorbed by the pulmo¬ 
nary blood in the alveoli. In contrast, only a very small 
portion (0.4%) of the inspired NO is absorbed by the 
airway. The remaining 44.5% of inspired NO is exhaled. 

In Fig. 10, Band C, the percentage of NO absorbed by 
the airway tissue is calculated under different condi¬ 
tions. In Fig. 10.6, Ci and flow rate change from 10 to 
100 ppm and from 50 to 1,000 ml/s, respectively. In Fig. 
10 C, simulations are performed twice, with and with- 


A 

Airway 

0.4% 



a Breatbhold 

~ j:rj ?5Q lav, (sec) 

v i (ml) 


Fig. JO. A: distribution of absorbed inhaled exogenous NO (Q - 10 
ppm) between alveoli and airway compartments for tidal breathing. 
B: effect of flow rate and inspiratory concentration on the percent NO 
absorbed by airway tissue (I^air)- C\ effect of inspiratory volume and 
breath-hold time on I % . ajr . 


A 



100 ppb 



Fig. II. Effect of reaction rate [or half-life on radial tissue 
concentration profile for 2 different gas phase concentrations; 100 
ppb {A) and 50 ppm { B ). 


out a 15-s breath hold, with different inspired volumes 
(Vi = 250-1,000 ml). The percentage of inspired NO 
absorbed by the airway tissue (I %ialr ) is inversely related 
to V E , Vj, and Vj and is a positive function of breath- 
hold time and Ci. Maximum absorption (8% of total NO 
inspired) in the airway tissue occurs when residence 
time in the airway compartment is maximized and 
contact with the alveolar compartment is minimized: Vi 
of 250 ml and 15-s breath hold. 

Reaction Rate Constant in Tissue 

To investigate the sensitivity of the exhalation profile 
on the reaction kinetics in the tissue, we varied E 
between 0.5 and 15 s, which corresponds to a range in k 
from 1.38 to 0.046 s -1 . Figure 11 demonstrates that Zj* 
values <4 s affect significantly the tissue concentration 
profile. Decreasing corresponds to a higher k or more 
rapid consumption of NO. The result is a reduced 
average concentration of NO in the tissue and an 
increased radial gradient at the airway wall. For low 
(Cj <K 418 ppb) inspired concentrations, the flux from 
the tissue to the airway is reduced. For high (Ci 418 
ppb) concentrations, the flux to the tissue from the 
airstream Increases. 

DISCUSSION 

Model Validation 

The model is able to simulate the different single- 
exhalation maneuvers (Fig. 3) extremely well. This 
finding becomes even more significant if one considers 
that only rough approximations of critical parameters 
were used (Dlno. S app , a i v , and S t , air ) that were acquired 
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by Independent means. No adjustment of model param¬ 
eters was made to optimize the model prediction of the 
experimental exhalation profile. The data are pre¬ 
sented for only a representative subject, inasmuch as 
our goal was only to simulate the fundamental features 
of the exhalation profile under dynamically changing 
conditions. The small differences between the experi¬ 
mental profiles and the simulations can be attributed 
to model simplifications and parameter estimation. 

The model assumes that the lung consists of two 
well-defined, separate regions: a rigid airway compart¬ 
ment and a well-mixed, expansile alveolar compart¬ 
ment. Such a division represents an idealization of the 
actual lungs. For example, the volume of the airways 
has been shown to change during breathing, there is a 
transitional region (respiratory bronchioles) between 
the airways and the alveolar region, and there may be 
axial concentration gradients within the alveolar re¬ 
gion (22), Although our model does not consider these 
features, the error introduced is minimal in relation to 
the model’s ability to describe the basic features of NO 
exchange. 

The model appears to fail to predict the exact value of 
the initial peak (phase I) in an inconsistent fashion. 
The peak is underestimated in Fig. 3, A, B, and D, and 
overestimated in the breath-hold maneuver (Fig. 3Q. 
Because our model does not include axial diffusion, one 
would expect the predicted peaks to be overestimated. 
Axial dispersion and mixing, especially at the interface 
between the alveolar and the airway regions, would 
create a broader peak width and hence a smaller peak 
height. The underestimation of the peak in Fig. 3, A, B, 
and D, is an artifact that arises from the fact that the 
model simulates a step change for the flow, from 
inspiration to expiration, while in reality a small period 
of decreasing Vj or zero flow might have occurred. In 
fact, a small breath-hold time (—0.5 s) between inspira¬ 
tion and expiration may represent more accurately a 
single-exhalation maneuver and wall increase the height 
of the peak in phase I. 

The model's simulation of phase III under the wide 
range of flow conditions is quite reasonable. Although 
we measured the average concentration in phase III 
experimentally over only a small range of flow rates, 
there is good agreement (Fig. 4). In addition, over a 
wider range of flow rates, the model is able to predict 
the same functional dependence reported by Silkoff et 
al. (24). 

The phase III slope under constant flow conditions 
and under dynamically changing flow conditions is 
similar to the experimental data (Fig. 5), but several 
important differences, exist. First, the intercept (phase 
III slope at constant Ve) is larger for the experimental 
data and importantly the model predicts a zero inter¬ 
cept after a J5-s breath hold. Second, the slope of 
b>iii ( No(r) vs. S nTiV _((j is steeper for the experimental 
data with and without a 15-s breath hold. Both of these 
discrepancies are likely due to simplifications in the 
model structure or parameter estimation. 

The model predicts a zero intercept after a breath 
hold simply because a period of 15 s is long enough for 


to reach its steady-state value (as discussed previ¬ 
ously). The structure of the model only allows for two 
mechanisms to create a phase III slope: 1) continuing 
gas exchange and 5) nonconstant absorption from the 
airway tissue. Continuing gas exchange is eliminated if 
C^., reaches Cai vss , and the model predicts that, under 
ambient or physiological conditions (C[ less than —50 
ppb), the absorption of NO from the airway tissue is 
essentially constant. Hence, the presence of a nonzero 
intercept experimentally after a 15-s breath hold is 
evidence that there are other mechanisms affecting gas 
exchange that cannot be completely neglected. For 
example, inhomogeneity in alveolar concentrations due 
to heterogeneous distributions of ventilation-to-volume 
ratios has been shown to contribute to the positive 
phase III slope of gases with similar physical character¬ 
istics (18, 19). In addition, axial or serial heterogeneity 
in the concentration of the gas in the respiratory region 
(respiratory bronchioles through the terminal alveoli) 
has been shown to contribute to the phase III slope of 
CO 2 , SF 6 . and He (22). The current model, of course, 
assumes a well-mixed alveolar compartment of uniform 
concentration that empties through a single path. The 
model also neglects axial diffusion, which may serve to 
increase Cai v above its steady-state value during a 
breath-hold maneuver. _ 

Regarding the slope of ,S liJ N o(') vs. S nl y E lO . a number 
of possible reasons can justify the small difference with 
the experimental data. Dl n0 has been demonstrated 
experimentally to be a positive function of lung volume 
(3). This contrasts with CO and is due to the relatively 
rapid kinetics of NO with hemoglobin. The rapid kinet¬ 
ics create a greater dependence of Dl no with surface 
area available for diffusion. The model assumes a 
constant Dl n0 throughout the exhalation. In addition, 
Dlno varies among subjects. The value of Dl n0 for 
subject 6 may be different from the mean value from 
Borland et al. (3) used in the simulations. The value of 
Dl n0 can affect the phase III slope by affecting the rate 
at which Cai v approaches C. a , vss [time constant is equal 
to V alv (i)/DL N0 ]. Inspiratory conditions affect the phase 
III slope (Fig. 7), and they were not controlled during 
the experiments. Finally, the phase III slope is curvilin¬ 
ear. Hence, the magnitude of 5) ir .NoM Is highly depen¬ 
dent on the interval of linear regression. 

A final point to consider related to the dependence of 
phase III on a changing flow rate is the relative 
magnitude of NO production in the airway and alveolar 
tissue. There is great intersubject variability in the 
slope and Intercept of the experimental data of F, NO vs. 
V E (Fig. 7 of Ref. 26); hence, there appears to be great 
variability in the relative ratio of NO production in the 
airway tissue to that in the alveolar tissue. For subject 
1, C a | VtSS is very small, yet the intercept is similar to that 
of subject 6. This is consistent with a larger relative 
rate of NO production in the airway tissueyhence, this 
subject demonstrates a steeper slope of .SnryicM vs. 
5m.v E (f) (Fig- 6 in Ref. 26). In the model simulations, 
values for S a[1[Kaly and S t , alr were chosen for a representa¬ 
tive subject but are only approximations. A small 
underestimation of S t ,ai r relative to S ap p,aiv can generate 
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a smaller slope in the ,S) nNO (f) vs. Sm.v r (f) diagram (i.e., 
less dependence of exhaled NO on the flow rate). 

Elimination Rate 

The positive dependence of Eno vs. V E has been 
experimentally reported (20, 24, 26). Silkof et al. (24), 
in agreement with our experimental findings (26), 
demonstrated a linear relationship between V E and E n q 
at the end of a single constant V E . In addition, we 
observed a similar dependence of E N0 in single exhala¬ 
tions with a dynamically changing V E . The model 
proposed by Hyde et al. (11) predicted an increase in 
E n0 with increasing ventilatory rate. With the assump¬ 
tion of constant lung volume, a steady-state approxima¬ 
tion, and inspiration of NO-free air, the}' demonstrated 
an increase in E N0 when the rate of ventilation was 
increased. Higher ventilatory rates resulted In a lower 
NO concentration in the lung, and thus the uptake of 
endogenously produced NO is shifted from absorption 
by the pulmonary blood to elimination by the exhaled 
breath. The above conclusion can satisfactorily describe 
a maneuver that has reached a periodic steady state 
(i.e., a steady-state rebreathing maneuver of NO-free 
air); however, this conclusion cannot be used to directly 
describe a single exhalation where the lung volume and 
concentration of NO are changing in time. 

Our model demonstrates that a nonzero alveolar 
concentration can satisfactorily describe the positive 
dependence of E N0 on V E . At the end of exhalation the 
alveolar concentration, which is independent of V E 
(expansile nature of alveoli, see Eq. 10), approaches a 
steady-state value; thus the positive dependence be¬ 
tween E no and V E becomes highly linear (Eq. 13a). 
However, there are at least two different possible 
explanations for the same phenomenon. First, Silkoff et 
al. (24) argued that at high flow rates the NO concentra¬ 
tion in the bulk gas will decrease, thus maximizing the 
gradient and the flux between the airway wall and the 
gas. Such an explanation is valid only if the transport of 
NO from the site of production to the bulk gas is limited 
by diffusion through the gas phase. Examination of the 
gas phase resistance along the airway (APPENDIX b) 
strongly suggests that, because of the very low solubil¬ 
ity of NO, resistance to diffusion in the gas phase is very 
small compared with diffusion in the tissue (6). Second, 
we suggested previously (26) that high flow rates would 
result in lower gas-phase and tissue-phase concentra¬ 
tions of NO and, hence, lower consumption by first- 
order chemical reactions and a larger flux of NO into 
the gas phase. However, the model demonstrates that 
first-order consumption of NO has only a minimal 
impact on the flux of NO at the airway wall. For 
example, if V B changes from 200 to 400 ml/s'in subject 6. 
E no increases by 75% (Fig. 6). However, <1% of the 
change in E no can be attributed to a decreased consump¬ 
tion in the tissue, even at the extreme of 0.5 s. 

The model’s ability to recreate the experimental 
relationship of p N0 to V E suggests that an experimental 
plot of E no vs. V e is an acceptable method of determin¬ 
ing the relative contribution of the alveolar tissue and 
airway tissue to exhaled NO. This has important 


implications for using exhaled NO as a noninvasive 
index of regional lung inflammation. For example, 
assessing airway inflammation in bronchial asthma 
might be accomplished by examining the magnitude of 
the peak in phase I. The main disadvantage of this 
technique is the difficulty in attaining a reproducible 
signal considering the short exhaled volume that is 
involved and the effect of ambient levels of NO on the 
magnitude of this peak. E N0 at end exhalation uses only 
the steady-state concentration of exhaled NO, which is 
independent of inspiratory conditions, yet when plotted 
against V E gives an accurate assessment of the relative 
production rates of NO in the airways and the alveoli. 
For example, we might anticipate an individual with 
active bronchial asthma to have a larger intercept 
(larger airway production) than normal (or compared 
with a previous test by the same subject before active 
inflammation) but the same slope. In contrast, an 
individual with pneumonia would have an increase in 
the slope but no change in the intercept. 

Exercise and hyperventilation have been demon¬ 
strated to increase Eno, but the mechanism is unknown 
(1, 17, 20, 21,25). Theories have included an increased 
production rate by the endothelial cells due to in¬ 
creased shear rate by the blood (17) and an increased 
Dl no (11). Our simple model suggests that the in¬ 
creased ventilatory rate associated with a nonzero C ah , 
can explain a large portion of the experimental observa¬ 
tions. The increased Dlno that is associated with 
exercise (4) cannot be used to explain the increased 
elimination, as suggested by Hyde et al. (11). In fact, an 
increased Dl no serves to decrease C^y (Eq. 12) and 
would thus tend to decrease E N0 . E N0 increases in 
response to exercise, because the increase in Dl n o does 
not lower Cai v enough to compensate for the increase in 
the ventilatory rate. 

Inspiratory Conditions 

Initial reports have suggested that ambient levels of 
NO do not affect the exhalation profile. Our model 
predicts that the NO exhalation profile, In particular 
the early part of phase III, is affected by not only the 
inspired concentration of NO, but also Vi and Vj. An 
interesting finding is the comparison of the peak values 
in phase I under the different inspiratory conditions. 
Inspiration of NO-free air obliterates the early peak of 
phase I in the absence of a breath hold. Similar 
experimental findings have been reported by Silkoff et 
al. (24), who reported no early peak when subjects 
inspired NO-free air. The absence of the phase I peak in 
these conditions does not indicate a lack of accumula¬ 
tion of NO in the airway, only that the accumulation 
leads to a concentration below C a i v . The model predic¬ 
tions related to inspiratory conditions provide a strong 
impetus and direction for further experimentation on 
the factors that control the NO exhalation profile. 

Exogenous NO 

Only a small amount (0.3-8%) of inspired exogenous 
NO at clinically relevant concentrations is absorbed by 
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the airway tissue. By use of a of 4 s, >90% of the NO 
that is absorbed by airway tissue is then absorbed by 
hemoglobin in bronchial circulation. This is evident by 
the relatively linear concentration profile of NO in the 
tissue layer at high Q. For a linear profile, the gradient 
(dC t /cLv) at the interface of the blood and tissue is the 
same as at the interface of the tissue and airstream; 
hence, the fluxes are the same: This model prediction 
depends on the value of the reaction rate constant, 
which may be underestimated at high tissue concentra¬ 
tions when reactions with 0 2 and thiols (2nd order in 
NO concentration) can no longer be neglected (2). In 
addition, the bronchial circulation is truly represented 
in vivo by discrete capillaries embedded with the tissue 
layer. Hence, some exogenous NO may diffuse through 
the tissue and reach important tissue layers such as the 
bronchial smooth muscle. Nonetheless, the current 
model provides an explanation for why exogenous NO 
has only a mild bronchodilator effect in humans (10): 
most of the exogenous NO that is inhaled is absorbed by 
the pulmonary circulation or is exhaled, and that which 
is absorbed by the airway tissue is primarily absorbed 
by the bronchial circulation before reaching the bron¬ 
chial smooth muscle. Furthermore, because a steady 
state in the tissue is reached in <1 s, the model predicts 
that the concentration of NO in the tissue depends only 
on Ci and not on rate of absorption. 

Conclusion 

A relatively simple and robust two-compartment 
model of NO has been developed that is capable of 
simulating many important features of NO exchange in 
the lungs. The model is able to simulate the experimen¬ 
tally observed dependence of exhaled NO concentration 
(and hence the dependence of the phase III slope and 
the elimination rate) on V E . This result stems from a 
nonnegligible contribution of NO in the exhaled breath 
from the expansile alveoli and the nonexpansile air¬ 
ways. The model predicts that continuing gas exchange 
in the alveoli can explain a large part of the phase III 
dynamics of NO. However, additional important param¬ 
eters that can influence the shape of the NO exhalation 
profile include preexpiratory breath hold, V[, Vj, and 
inspiratory NO concentration. The effect of inspiratory 
conditions on the exhalation profile has not been well 
studied experimentally but may play a critical role in 
the the interpretation of the NO exhalation profile. 
However, the model also predicts that the alveolar 
concentration approaches its steady-state value if exha¬ 
lation is >8 s long. Thus the expiratory concentration is 
independent of inspiratory conditions if duration of 
exhalation is greater than ~8 s. Finally, the model 
suggests that the relationship between En 0 at end- 
exhalation (exhalation time > 8 s) and V E may be a 
simple, effective, and reproducible technique for deter¬ 
mining the relative contribution of the airways and 
alveoli to exhaled NO. 

APPENDIX A 

The unsteady diffusion equation for the tissue compart¬ 
ment {Eq. 1) is solved for the simplified case of a steady state 


(Eq, 3). To justify this assumption, the unsteady equation 
must be solved and analyzed. If the gas phase concentration 
suddenly increases from 0 to a constant value Ci (i.e., Inspired 
concentration), the concentration profile of NO in the tissue 
will be given from the solution of Eq. 1 with the following 
boundary conditions 


C t (r, 0) = 0 (Ala) 

C t (f. L Ciair ) = C w = CA.-g (Alb) 

The initial condition is derived from Eq. 3 for C w = 0 


c,<p, 4 -«*%***-<*%* 


-x/fS . 


Aair 

T" 


(Ale) 


where 


— 2 -j- 

“i = Pw?. e-L.we ' “2=-“i +1 


In dimensionless form, Eq. 1 can be written as follows 



— - p 2 © = — 
dz? d't 

(A2a) 

with 

0(0, z) - cue* 1 * - cL 2 e~ vz 

(A2b) 


0(r, 0) = -1 

(A2c ) 

where 

Cw air 

©(t, 1) = -- = £ 

S fk 

(A2d) 

q A - S t , a | r /& 

X 

t 

F = hj 



AU'A ’ 


By use of the method of separation of variables, Eq. A2 can be 
solved to give the following solution 


x. -< iJin 2 2 2 

0 = Zu ~ 7 —"T - ; r - & T ) sin (rm$ 

n-l \1 Z + Fr r i\ 


ZrtiiL, 


- (e-' ,V - - 


jSTl JJL 2 + /J 2 'TT 2 


d*~) sin (mr( 1 — z)l ( A3) 


+ 2 


2nrr[(— 1)" - 1] „ 

_——_ 1 3“ /TIT T 

p 2 + it-rE 


sin (rmz) 


In Fig. 12 the solution of the dimensionless concentration (©) 
is plotted as a function of the dimensionless radial distance (z) 
for different values of the dimensionless time (t). The profiles 
computed from Eq. A3 for p z = 0.524 and l = -0.97 
correspond to a tissue thickness of 100 pm and a step change 
in the gas concentration of NO from 0 to —60 ppb, respec¬ 
tively. One can see that atr - 0.2 (f = 0.6 s) the concentration 
of NO at any given point in the tissue is within 10% of tile 
steady-state value. As a result, under normal conditions (i.e., 
inspiration of room air of small NO concentration and expira¬ 
tion time >2 s) the error introduced by the steady-state 
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From Eq. Cl 

V alv (i) = v aW (0) - v Ep0 f - sm (C2) 

Then, the solution to Eq. 10 will be, for case I. for A = V| 0 + 
4sV a |„{0) > 0 (always true when the flow rate is increasing) 


CalvM ” 


c 

'"'□pp.alv 

Dl N o 
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DU, 


-^app.^lv 


Dl, 


NQ 


^,-[2DLN6Tanh “ 1 [(V^o+ZsfV t/AJ/^A} ( C3s) 

CDlno 


c alv (0) - 


-'app.alv 


Dl, 


NO 


Fig, 12. Effect of a step change in gas-phase concentration on radial 
tissue concentration profile from 0 to 60 ppb. 


assumption is negligible. However, for a fast breathing maneu¬ 
ver of a concentrated gas mixture of NO, such an assumption 
is not valid. 


where 


h. 2 ~ 


-Ve,6 ± >/A 

2s 


for case II, For 
A = V| 0 + 4<rV alv (0) < 0 


APPENDIX B 

The radial gas phase resistance for the diffusion of NO is 
negligible and can be proven using the following argument: If 
there is a significant gas phase resistance in the transport of 
NO between the tissue and the bulk gas, then Eq. 4 can be 
rewritten as 


—^tig.sur ® 6C W f'sir) 

and solving for Cw 


(Bl) 


decreasing flow rates with |sj > 


V|,o 

4V aIV (0). 


C 3 l v (I) 


C 

DL no 


+ 


C a ,v(0) 


Cl _ 

a -i2(DL N o)Tan- , [(-V E ,o-2iDACil/v'WI 

dl no J 


{C3b) 


_ & + -KsC atr 

,/J ~ b + KJK c , & 


(B2> 


where K s is the gas phase mass transfer coefficient for NO, \ t;g 
is the partition coefficient between the liquid and the gas 
phase, C air is the the concentration of NO in the bulk gas, and 
a and b are constants from Eq. 4. 

The Sherwood number (Sh = K a d alr /D SS D in the airway 
tree, although dependent on the flow rate and the axial 
position, has values usually between 6 and 10 (6). The 
diffusivity of NO in the gas phase (D gas ) is 0.23 cm 2 /s, and the 
diameter of the airway tree (d alr ) is between 2 cm and 5 mm. 
Thus, even under the extreme conditions, K s has a minimum 
value of 0.69 cm/s. The values for the constants a and b, as 
estimated for the typical subject of our simulations, are 0.067 
ppb-cin -1 -s -1 and 0.0039 cm/s, respectively. Because (KJK^ » b 
and (atIQ <s C air , Eq. B2 can be written as 


APPENDIX D 

The mass balance in the upper airway compartment can 
also be solved for nonconstant flow rates. For the simple case 
of a linearly changing V E and a linear function of V, Eq. 

6 becomes 


* v c/-i |r 

-CVb.0 + St) — + (a - ^«ir) (% + beV) = — (Dla) 
with the general initial and boundary conditions 

c alr . ( „ 0) = m ( Dib ) 

C a lr.{V-0) “ {Elcj 

In dimensionless form, Eq. Dl can be written as 


~ \:g C alr C B 3) 

Hence, the NO exchange between the tissue and the airway is 
not gas-phase diffusion limited. 

APPENDIX C 

Solution for the mass balance of NO in die alveolar 
compartment (Eq. 10) can be_ easily derived for a simple case 
of linearly time-dependent V E . For an initial flow rate and 
slope (Ve.o and s, respectively) 

V E « = V Ei0 + st (Cl) 


da da 

— + (2 + 2t) — = ~(c + dz)£ 
At flz 


(D2a) 
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where 


C a[r - a!b\,. 

Vee^-s 
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fiV E . o 

c — - a E , d = 


V E . 0 /s V| 0 /(2s) 

^l.o^e 


2.S 2 


Eq. D2 can be solved using the method of characteristics to 
obtain region 1 for (t + l) 2 - I < z 

/! — [Vl 0 /(2s)] nl — al{bk,. J - 

5(f, ^ -- - (D3a) 


4 


region II for (t + l) 2 - I > z 

gjV E u/s{y!u + 1 - 1 )] - a/{b\ l:g ) 

V E . 0 /s {D3b) 

x gd v^T- i ~d + ^ < m. 

with 

u = (f + I) 2 — 1 — z 

A = L M\' U + 1 — l) 2 — u(>/i/ 4- 1 — 1) + tu — V^t(t + z) 
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